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Abstract 

A new method to study the long-range correlations in multiparticle 
production is developped. It is proposed to measure the joint factorial 
moments or cumulants of multiplicity distribution in several (more 
than two) bins. It is shown that this step dramatically increases the 
discriminative power of data. 

1 Introduction 

It is now widely recognized that long-range correlations (LRC) in rapidity 
give information about the early stages of the collision. Indeed, such corre- 
lations cannot appear at late stages in the evolution of the produced system 
when longitudinal expansion separated the particles by large distances. Just 
after the collision, however, the system is small enough for the correlations 
to extend through the whole system. 
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A special case of LRC are forward-backward correlations where one com- 
pares particle distributions in two intervals located symmetrically in the for- 
ward and backward hemispheres. They were extensively studied since early 
times of high-energy physics PQ. 

One of the interesting issues in particle production is the question if the 
produced particles "remember" the colliding projectiles, their energies, mo- 
menta and quantum numbers. Obviously the answer depends on the kine- 
matic region we are considering. Close to the fragmentation region, the 
influence of the projectile on the produced particle spectrum is naturally ex- 
pected. In the central rapidity region, far from the projectile fragmentation, 
the question remains open. On the theoretical side there is no consensus and 
various models give different answers. 

An excellent review of models can be found in |T] (see also, j2]), therefore 
here we only quote some examples. With respect to the question of the 
number and structure of particle sources, they may be divided into three 
categories. In the first one, originating from the famous Landau and Feynman 
papers [21 H] , particles produced in the central rapidity region are decoupled 
from the projectiles. Thus the source of particles is symmetric with respect 
to y = 0. In the second class, like the wounded nucleon model [5], particles 
are produced by quasi-independent emisssion from the two colliding objects. 
In this case particles in the central region come from two sources, naturally 
asymmetric ones [HI E] . There is of course also a third class which combines 
the two pictures, a typical example being the dual-parton model [5]. 

These various mechanisms can be tested (and verified) by studying the 
forward-backward correlations. The essential point is that correlations for 
one symmetric source are generally much stronger than those induced by 
two asymmetric ones (21 IS] . Following this general idea we recently proposed 
a systematic method of investigation of the forward-backward correlations 
in symmetric hadronic and heavy ion collisions [2] . It was shown that such 
investigations allow to verify how many indepedent sources of particles con- 
tribute to the observed distributions. 

In the present paper we generalize these results in two respects: 

(i) we abandon the requirement of symmetry and consider the general 
case of asymmetric processes and thus also asymmetric sources; 

(ii) We suggest to measure and compare particle distributions in more 
than two intervals^. 

1 Measurements in three intervals were considered in [10] and [llj . 
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This generalization of the problem allows to undertake a general discus- 
sion of LRC and thus extends its application to other processes, like, e.g. 
lepton-nucleon, hadron-nucleus and asymmetric nucleus-nucleus collisions. 

We consider measurements of multiplicity moments in B intervals. Fol- 
lowing our previous paper [2J, we assume that particles are produced by 
indepedent sources and that the population of particles in B bins from one 
source is random, i.e. it is described by the multinomial distribution!!. 

We then evaluate the number of measurable (factorial) moments of the 
distribution and compare it with the number of parameters in the system. 
This allows to draw our main conclusion: the discriminating power of the 
method increases dramatically with increasing number of intervals in which 
the measurements are performed. 

In the next section we present the mathematical structure of our ap- 
proach. In Section 3 the number of possible measurements and number of 
independent parameters are evaluated for the general case of particle produc- 
tion from independent sources and a measurement in B intervals. An impor- 
tant special case, when the number of independent moments is reduced by 
symmetry, is discussed in the Appendix. In Section 4 the explicit formulae 
giving the factorial cumulants for two models with fixed numbers of sources 
are derived and commented upon. The general formulae and two examples of 
models with fluctuating numbers of sources are discussed in Section 5. The 
summary of the results can be found in the last section. 



2 Formulation of the problem 

Following the assumptions explained in Introduction, we write the generating 
function for the particle distributions in the B bins in the form 

$(*!,..., Zfl) = (j[4>T i {pHZ 1 + ...+pBiZ B )), (1) 

where 0, is the generating function for the i-th source, pki is the probability 
that the i-th source sends particles into the k-th bin and Wi is the number of 
sources of type i. The angular brackets denote averaging over the multiplici- 
ties Wi (they can be omitted if the W{ do not fluctuate). Since the generating 

2 This assumption which may be understood as the definition of a "source" is accepted 
in practically all published discussions of the forward-backward correlations, see e.g. |12j . 
For an extensive list of references, see pQ. 
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functions 0j do not have to be all different, one can assume without loss of 
generality that each of the numbers Wi can only take the value zero or one 
(Section 5.1), but sometimes it is more convenient to assume that Wi can be 
any nonnegative integer (Section 5.2). The measurable (factorial) moments 
are given by 

p. -/lT n r \ - d r ®{z u ...,z B ) 

where rij is the number of particles in bin j. Here and henceforth all deriva- 
tives are taken at z — Z\ — ...zb = 1. Note that 

dT ^ Z ) = F (r) (1) 

is the r-th factorial moment of the distribution of the total number of particles 
sent by source n to all bins. 

If the numbers and nature of the sources do not fluctuate, it is advanta- 
geous to introduce the cumulants 



d r log $(z u . . . , z B 
dzl 1 . . . dz*i 

which, as is easily derived from (TI]), can be expressed as 

N 

fh...i B = E Pin ■ ■ -Pb ti\ (5) 

n=l 

where 

f (r) = dfog(j>n(z) . . 

Jn ^ \») 

are cumulants of the distribution produced by the n-th source. 
Let us also note here that using (JHJ) and the identity 

E ~ r T~~ r ~\P^x ■ • -Pb 9 = (Pi + -+PbY = 1, (7) 
3v--3b\ 

one finds the very useful relation 

r! N r 

Jv--JbI n= i 
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3 Counting of parameters 



Consider a general situation of N groups of independent sources, all sources 
in one group being identical, and B bins. No symmetry relations among 
groups are assumed. 

We first evaluate the number of moments which can be measured. To 
this end we observe that each moment has B indices: Fi^ 2 ___i B . Define the 
rank r of the moment as 



t = h + 12 + ■■■ + is- (9) 

The number of moments at given r and B, m(r,B), is the solution of the 
well-known combinatorial problem: in how many ways can one distribute r 
identical objects among B boxes: 



(r + £-!)! r — (B + r max )\ 

m{r > B) = ^U^W - g^'^-Bi^r- 1 (10) 

The next thing we want to know is the number of parameters in the 
model. First, there are N(B — 1) independent probabilities. In addition we 
need also, for each kind source, the derivatives of order up to r rnax of the 
multiplicity generating function [(p n (z)} Wn : 

^S r) = ^[C*(*)]- (ii) 

They are polynomials in the random variable w n . 

The expressions for the measurable moments of order r contain the aver- 
ages 

N 

(F^...F^), £r„ = r. (12) 

n=l 

When the multiplicity distribution for sources is not known, each of these 
averages is an independent parameter. Using the same combinatorial for- 
mulas as before, we thus find that the number of independent parameters 
is 

P(B,N,r max ) = N(B - 1) + - 1. (13) 

'"max-'* • 
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When the distribution of numbers of sources W(wi, . . . ,wn) is known, 
all averages (1121) are determined in terms of F^™' 1 and therefore the number 
of independent parameters is 

P{B,N,r max ) = N{B + r max -l). (14) 

Thus we finally obtain for the number of parameter-independent con- 
straints between the measurable quantities 

C(N; B; r max ) = ^ + r — ) ! _ i _ P(jB; Nj w ) (15) 

where P(B, N,r max ) is given by f fT3|) or ffT4"|) . 

To obtain tests, we demand that C > 1. It is clear that for any N and 
r m ax > 2 one can always find B such that this condition is satisfied. 

For practical reasons, one has to keep r max rather small, say 2 or 3. In 
Table 1 we give the minimal number of bins necessary to obtain parameter- 
independent constraints. 

Table 1 
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To illustrate possible applications of this general discussion, we present in 
the next two sections four examples of specific models of particle production 
which can be tested in this way. 

4 Fixed number of sources 

For a fixed number of sources the measurable cumulants are given by ([5]). 
Below we give two specific examples. 
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4.1 Landau model: one source 

In the Landau model there is just one source of particles, resulting from 
hydrodynamic expansion of the remnant of the two projectiles (just after 
collision the remnant is concentrated at y cm = 0). For one source (N = 1) 
already at r max = 1 the number of measurable moments is equal to the num- 
ber of parameters. Therefore, it is possible to determine all the probabilities 
Pj from the moments (cumulants) of rank one. For each r > 1 there is one 
more parameter, f( r > = d r [log (p(z)]/dz r . Using fl5]), this parameter can be 
evaluated from any measured moment of rank r. Indeed, for one source we 
simply have 

f£L, jB pf..*g7 (p) (i6) 

where the subscript denoting the source was dropped. Since all probabilities 
are already determined from the moments of rank one, this formula allows 
to evaluate and thus all other measurable cumulants of rank r. 

4.2 Deep inelastic scattering: two sources 

In deep inelastic scattering there are at least two different sources: the proton 
and photon remnants. It is thus interesting to investigate if these two sources 
are sufficient to describe the data. In this section we show that the hypothesis 
of two sources gives indeed strong constraints on particle correlations. 

Following the argument of Section 3, we consider B bins located anywhere 
along the direction of the incident photon. We thus have 2(B + r max — 1) 
parameters (the probabilities pj 7 , pjp and the cumulants ffl, fp ). Let us 

(r) 

denote by fj the measurable cumulant of order r of the distribution of 
particles in the bin j. We show below how, using the measured cumulants 
fj V ^ for r < 2, one can determine all the probabilities. 

Since the sum rule ([H]) allows to determine the sum of the cumulants 
/+ = + fp^ for any r, we are left with with r free parameters = 

ffl — f P which should be sufficient to predict the correct values of the other 
measured cumulants. Instead of the parameters pj 7 pjp it is more convenient 
to use 

Pj±= Pn±PjP- ( 17 ) 
As already mentioned, for each r the parameter can be obtained 
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directly from formula (jSJ). For r = 1 and r = 2 we have 



B B B 

f^=Y / <n l >; /j 2) =Y i <n i (n i -l)>+2Y l <n i n j > (18) 

8=1 8=1 i>J 

where is the number of particles observed in bin i. 

Let us consider first the cumulants of order one (they coincide with the 
moments of order one, i.e. average multiplicities). From formula (J5J) one gets 

which together with the sum rule for eliminates B parameters. 
Let us consider now the cumulants /j . From (jSJ) we have 

4/f = P?_/f + 2 Mj+ /i 2) + (20) 

Using (|T9|) to eliminate pj+ we get a quadratic equation for The two 
solution of this equation depend on the parameters and Thus we 
get 2 B possible sets pi_ . . .ps-- Hopefully most of them can be eliminated 
by the obvious requirement that each pj_ must be real and that the following 
constraints must be satisfied. 

B 

\pj-\ < p j+ , \Pj-\ < i, J2pj- = °- ( 21 ) 

3=1 

Thus, if the model is consistent with data, i.e. solutions exist, all the proba- 
bilities are determined, though some ambiguities may be left. 

For r max = 2 we have, in addition, \B(B — 1) cumulants of the type 

/110...0 which should be fitted with two parameters and /£ . Increasing 
fmax by one, introduces two new parameters /+ and /i 3 \ The former, 
however, is fixed by the sum rule ([8]) so that there are m(3, B) new cumulants, 
constrained by the sum rule which has already been used, to be fitted with 
one free parameter. 



5 Fluctuating number of sources 

When the number and nature of sources fluctuate, the discussion of LRC 
becomes rather involved. The reason is that the formulae expressing the 
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measurable moments in terms of the parameters of the model become com- 
plicated, as can be seen later in this section. In most models, however, the 
sources are not entirely arbitrary and thus these relations can be simplified. 

We start with the general formulae for arbitrary number and nature of 
sources and then discuss two examples, suggested respectively by the dual 
parton model and by the wounded constituent model. 

5.1 General formulae 

Let us consider the generating function (JTj) with each Wi equal zero or one. 
Then 

Wi (wi - 1) = 0. (22) 
This greatly simplifies the differentiations. In fact 

^4>r{z)=w<^M*)=Vi*i r) . (23) 

Below we give the formulae for the measurable factorial moments of rank 
1, 2 and 3. They are written assuming that only the first (for r = 1), the first 
two (for r = 2), or the first three (for r = 3) bin indices are non- vanishing. 
Analogous formulae are of course valid for any other selection of bins, pairs 
of bins and triplets of bins. 



N 

F w ... = "£(^i)PuF^\ 

i=i 

N N 

F 110 ... = J2( w i)PuP* F i 2) + J2(^j)PuP2jF^F^ (24) 
i=i i^j 

N N 

Fmo... = J2( W i)PliP2iP3i F i + J2 (WiWj)PliP2iP3jF[ 2) F!> ^ + 
i=l i^j 
N 

E (mw jWk ) PliP2jP3k F^F^Ft'\ 

(r) 

where F{ is the r-th factorial moment of the distribution of particle from 
source i (c.f. (j3J)). When some indices coincide, it is enough to change 
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correspondingly the bin indices of the probabilities p. For instance, 

N N 

F20... = J2( w i)P 2 u F i 2) + E<Wi)PHPi^ (1)i f } - (25) 

i=l ijtj 

5.2 Dual parton model 

For a general nucleus-nucleus collision we have a certain number Nl of iden- 
tical sources moving left, a number Nr of identical sources moving right and 
Nq identical symmetric sources. These numbers fluctuate from event to event 
and their distribution depends also on the centrality of the collision. The left 
and right moving sources are mirror images of each other with respect to cm 
rapidity. 

We consider the case where the bins are also selected to be symmetric 
with respect to y cm = 0. Then if (fi a (p± a Zi + . . . + PbuZb), where a stands 
for asymmetric, is the generating function for the multiplicity distributions 
in the bins 1, . . . , B of the particles originating from a left moving source, 
then 4> a (PBaZi + • • • +Pio,zb) is the corresponding generating function for the 
particles originating from a right moving source. 

Let us denote by w^, Wr, Wq the numbers of left moving, right moving 
and central sources. In j2] we discussed mostly the case of two bins and 
fixed wl = wr and wc- Here we assume an arbitrary number of bins and 
a general joint probability distribution W(wl, wr, wq) which, however, can 
be evaluated, e.g. by the Glauber method (the result will, naturally, depend 
on the model adopted for particle production). Then the overall generating 
function for the multiplicity distributions in the B bins is 



...,z B ) = W{w L ,w R ,w c ) 

[4>a(PlaZl + ••• +PBaZEi)] WL [4>a{PBaZl + ••• + Pl a Z B )] WR 

[MpicZi + -+Pbcz b )] Wc ■ (26) 

We will denote the probabilities by p, L A where A = L,R, C . Although the 
probabilities piR can be expressed by the probabilities pn, this redundancy 
in the notation makes the following formulae much shorter. Similarly, the 
derivatives ( ITT]) are denoted by Fj[\ Using this notation, the explicit ex- 
pressions for the measurable factorial moments, obtained by differentiation 
of flU, read 
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.4 




F%l=Y.PiAPjAP k A(FT) + 



A 



+3 E [PiAPjBPkB + PjAPkBPiB + PkAPiBPjB}{F^ F^} + 



A^B 



+ E PiAPjBPkciF^FUpV). 



(27) 



A±B±C±A 



The parameters of the model are the 2(5 — 1) probabilities and the av- 
erages ffl2l) . Their number is given by ffl3|) or (}Tlj) with N = 2. 

Since formulae (1271) are rather complicated, it seems that in absence of 
other constraints, the best way to proceed is to try to fit them by minimizing 
the x 2 - If th e fit works, the resulting values of the probabilities and of the 
factorial moments give information about the properties of the sources. 

For pp scattering wl = w R = 1 and thus the relations are much simpler. 
As they may be easily obtained from (I27p . we discuss here only the number 
of parameter-indepedent constraints. The number of possible measurements 
is given the the Appendix. The number of parameters is I? + -B/2 — 2 + 2r max 
for B even and B + (B + l)/2 — 2 + 2r max for B odd. One can see that 
for r max = 3 there are already 3 constraints for B = 3 and 8 constraints for 
B = 4. If one wants to restrict the measurements to r max = 2, it is necessary 
to measure distributions in at least 5 bins. Then one obtains 2 constraints. 

5.3 Wounded constituent model 

In the wounded constituent model, particles are emitted by the wounded 
constituents moving left or right, thus there are no central sources. The 
relevant formulae can be obtained from ( |27|) by putting F^ = 0. 
They can be written in the form: 



,(2) 



FP=P lL {F^)+p lR {F^)- 

PiLPjL(Fi 2) )+p lR p jR (FR) + (PiLPjR + PiRPj^iF^F^); 

F$ = {p l LP j LPkL(F^ ) ) + (L^R)} + 
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+ {3[PiLPjRPkR + PjLPkRPiR + PkLPiRPjR] (F^ F R 2) ) + (L 4+ R)j . (28) 

The consequences of these formulae are different for symmetric (e.g. Au — 
Au collision) and asymmetric (e.g. d — Au collision) processes. 

When the probability distribution W{wl,Wr) is known, the number of 
parameters at a given r max is B + r max — 1. For asymmetric processes the 
number of possible measurements is given by (fit)]) and for symmetric pro- 
cesses the relevant formulae are given in the Appendix. Already for B = 3 
and r max = 2 one obtains 2 parameter-indepedent constraints for symmetric 
and 5 for asymmetric processes. 

When W(wl,wr) is not known, the various averages of the moments 
(Fl Fr ^), (s = 0, ...,r), have to be fitted from data at every r < r max . At 
given r the number of independent averages is r + 1 for the asymmetric case 
and for the symmetric case it is r/2 + 1 for r even and (r + l)/2 for r odd. 
It is remarkable that already at B = 3 parameter-indepedent constraints 
exist. For symmetric processes one obtains 1 constraint for r max = 2 and 5 
constaints for r max = 3. For asymmetric processes the correponding numbers 
are 2 and 8. 

6 Summary 

Extending the ideas formulated in [2] (see also [9]), a new method to study 
the long-range correlations in particle production is developed. The new 
proposition is to measure the factorial moments and/or cumulants in several 
bins, as opposed to previous studies which were mostly restricted to just two 
bins (see, however, footnote 1). It was shown that increasing the number of 
bins magnifies dramatically the possibility of discriminating between various 
models of particle production. 

The discriminative power of the method was analyzed in the most general 
way. Apart from this general treatment, four specific (and popular) models 
of particle production were discussed. It was shown that the suggested mea- 
surements provide strong constraints on all of them. 

The method seems rather general and flexible. It can be applied to sym- 
metric, as well as to asymmetric processes. It can be used to study distribu- 
tions in various kinematic variables (e.g. rapidity and transverse momentum 
[2]). Finally, it does not require full acceptance of the detector. 
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We conclude that studies of long-range correlations in multiparticle pro- 
duction may become a powerful instrument in investigations of particle pro- 
duction mechanisms at high energy. 

7 Appendix 

In this appendix we calculate the number of moments of order r for the 
reflection symmetric case, i.e. when 

Fil...i B = ^iB—ii- (^9) 

We will denote the number of these moments by ms(r, B). Let us call sym- 
metric the moments for which the ordered sets {ii-.-is} an d {%b ■ ■ -h} 
coincide. The number of such moments will be denoted S(r,B). The con- 
straint (129]) does not affect the number of symmetric moments, but reduces 
the number of the other independent moments by a factor of two. Thus 

m s (r, B) = i[m(r, B) + S(r, B)], (30) 

where m(r, B) is given by formula (fTUj) and S(r, B) remains to be calculated. 

Let us begin by the case when B = 2K + 1, where K is an integer. 
Then the generic form of a symmetric moment is i ? i 1 ,...,i K , n ,i iC ,...,ii- Therefore, 
the number of such moments can be calculated as follows. Include all the 
nonnegative integers p such that n = r — 2p > 0. Notice that for each p there 
are m(p, K) moments, thus in this case 

%^) = E( P+ p _1 )- (31) 

The case B = 2K reduces to the previous one with the constraint that n = 0. 
For r = 21 + 1, where I i an integer, there are no solutions for p, therefore 

S{r,B) = Q, (32) 
while for r = 21 the only solution is p = I, so that 

S(r,B)=( I + K I (33) 
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